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Quantum  mechanical  calculations  are  used  to  predict  gas,  liquid,  and  solid  heats  of 
formation  of  energetic  molecules.  A  simple  atom-equivalent  method  converts  quantum 
mechanical  energies  of  molecules  and  their  atomic  constituents  to  gas-phase  heats  of 
formation  of  energetic  materials.  Functional  relationships  between  heats  of  vaporization 
and  sublimation  and  properties  associated  with  quantum,  mechanically  determined 
electrostatic  potentials  of  isolated  molecules  are  established.  These  are  used  with  the 
gas-phase  heats  of  formation  to  predict  condensed-phase  heats  of  formation.  The 
calculated  gas-phase  heats  of  formation  have  a  root  mean  square  (rms)  deviation  of 
3.1  kcahmol  and  a  maximum  deviation  of  7.3  kcal/mol  from  35  experimental  values. 
The  rms  and  maxirmim  deviation  of  predicted  heats  of  vaporization  from  27  experimental 
values  are  1.7  and  6.1  kcal/mol,  respectively.  The  rms  and  maximum  deviations  of 
predicted  heats  of  sublimation  from  36  experimental  values  are  3.6  and  12.4  kcal/mol, 
respectively.  The  rms  and  maximum  deviations  of  predictions  of  liquid  heats  of 
formation  from  41  measured  values  (corresponding  to  24  molecules)  are  3.3  and 
9.3  kcal/mol,  respectively.  Similarly,  the  rms  and  maximum  deviations  of  predictions  of 
solid  heats  of  formation  from  75  measured  values  (corresponding  to  44  molecules)  are 
9.0  and  35.4  kcal/mol,  respectively. 
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INTRODUCTION 

The  availability  of  propellants  with  significantly 
improved  impetus  and  controlled  burning  rate 
properties  is  essential  to  the  successful  develop¬ 
ment  and  deployment  of  new  high-performance 
tank,  artillery,  electrothermal- chemical,  and 
Naval  bombardment  systems.  Design  concepts 
for  the  formulation  of  these  advanced  propel¬ 
lants  is  a  crucial  emerging  technology  in  the 
United  States  Department  of  Defense,  where 
limited  funds  dictate  optimization  of  time  and 
resources.  In  previous  times,  formulators  would 
develop  propellant  mixes  in  a  series  of  trial-and- 
error  formulations,  and  heats  of  formation  of 
the  ingredients  (if  unknown)  were  measured  for 
use  in  calculations  to  predict  the  impetus  of  the 
propellant  under  gun  firing  conditions.  The 
determination  of  impetus  is  one  of  the  first  steps 
in  the  screening  process  of  propellant  formula¬ 
tion.  Candidate  propellants  that  show  enhanced 
performance  receive  additional  study;  poor  per¬ 
formers  are  eliminated  from  further  consider¬ 
ation.  Waste  emanating  from  measurement  of 
heats  of  formation  of  unsuitable  candidates  can 
be  reduced  through  development  and  use  of 


’Corresponding  author:  Dr.  B.  M.  Rice,  U.S.  Army  Re¬ 
search  Laboratory,  AMSRL-WM-BD,  Aberdeen  Proving 
Ground,  MD  21005-5066.  E-mail:  betsyr@arl.mil 

COMBUSTION  AND  FLAME  118:445-458  (1999) 

©  1999  by  The  Combustion  Institute 
Published  by  Elsevier  Science  Inc. 


theoretical  capabilities  for  prediction  of  this 
property.  The  research  presented  here  describes 
such  a  development  of  theoretical  capabilities 
to  aid  in  the  screening  process  of  propellant 
formulation.  In  this  work,  we  present  computa¬ 
tional  tools  that  convert  quantum  mechanical 
calculations  of  energetic  materials  to  heats  of 
formation  in  the  gas,  liquid,  and  solid  phases. 

There  are  a  variety  of  methods  to  predict 
gas-phase  heats  of  formation  from  quantum 
mechanical  information.  One  method  uses 
known  heats  of  formation  of  isolated  atoms  and 
calculated  atomization  energies  {Dq)  to  predict 
gas-phase  heats  of  formation  of  molecules  [1]. 

L^iA^y,  OiO  =  x^^{A,  OK) 

+  yA/{^{B,0K)-2Do 

(1) 

This  method  was  found  to  reasonably  predict 
the  heats  of  formation  for  a  variety  of  organic 
and  inorganic  molecules  with  the  best  predic¬ 
tions  corresponding  to  those  using  the  G2  level 
of  theory  [2, 3].  While  the  degree  of  accuracy  of 
the  predictions  using  this  level  of  theory  has 
been  impressive,  the  calculations  require  com¬ 
putationally  expensive  electron  correlation 
treatments  which  might  be  prohibitive  for  sys¬ 
tems  containing  a  large  number  of  atoms  or 
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where  computational  resources  are  limited.  The 
study  also  assessed  the  performance  of  the  more 
computationally  tractable  density  functional 
theory  (DFT)  [4-7]  in  calculating  heats  of  for¬ 
mation  [1,  8].  The  DFT  predictions  (using  the 
B3LYP  density  functional)  [9],  while  not  as 
good  as  the  G2  predictions,  are  reasonable.  A 
similar  procedure  for  predicting  heats  of  forma¬ 
tion  using  a  modest  level  of  quantum  mechani¬ 
cal  theory  has  been  shown  to  be  accurate  to  a 
few  kcal/mol  [10].  The  method,  known  as  BAC- 
MP4,  corrects  for  errors  in  the  level  of  theory 
through  empirical  bond-additivity  corrections. 
Melius  applied  this  method  to  90  molecular 
species,  and  reports  an  average  deviation  from 
experiment  of  1.3  kcal/mol  [10]. 

Another  method  of  predicting  gas-phase 
heats  of  formation  is  based  on  Hess’  Law  [11] 
and  uses  a  combination  of  quantum  mechanical 
and  experimental  information.  Hess’  Law  states 
that  the  standard  reaction  enthalpy  is  expressed 
as: 

AF/keaction  =  (ProduCts) 

—  A//^  (Reactants)  (2) 

The  standard  heat  of  formation  of  a  single 
component  of  a  reaction  (either  product  or 
reactant)  can  be  determined  using  the  reaction 
enthalpy,  which  can  be  obtained  from  quantum 
mechanical  calculations,  and  reliable  values  of 
heats  of  formation  of  the  remaining  products 
and  reactants.  The  accuracy  of  the  prediction 
using  Eq.  2  is  increased  if  the  reaction  is 
isodesmic  [12].  Isodesmic  reactions,  in  which 
numbers  of  electron  pairs  and  chemical  bond 
types  are  conserved  in  the  reaction,  allow  for  a 
cancellation  of  errors  inherent  in  the  approxi¬ 
mate  treatment  of  electron  correlation  in  the 
solutions  to  the  quantum  mechanical  equations. 
To  illustrate,  we  predicted  gas-phase  heats  of 
formation  of  l,3,5-trinitro-l,3,5-s-triazine 
(RDX)  (experimental  value  =  45.8  kcal/mol 
[13])  using  Eq.  2.  We  calculated  the  heat  of 
reaction  of 

3CH2NNO2  ^  RDX 

using  DFT  at  the  B3LYP/6-31G*  level  [9,  10, 
14]  and  combined  this  with  a  reported  heat  of 
formation  of  the  reactant,  CH2NNO2  (33.6  kcal/ 
mol)  [15].  The  predicted  heat  of  formation  for 


RDX  is  52.8  kcal/mol,  ~7  kcal/mol  greater  than 
the  experimental  value.  This  method  of  calcula¬ 
tion,  while  giving  a  reasonable  prediction,  re¬ 
quires  reliable  values  of  heats  of  formation  of 
other  components  of  the  reaction. 

Atom  equivalent  schemes  are  used  to  convert 
quantum  mechanical  energies  of  formation  of 
atoms  to  heats  of  formation  for  various  classes 
of  molecules.  The  gas-phase  heat  of  formation 
using  atom  equivalents  is  represented  as: 

^Hi  =  E,  -  S  njej  (3) 

where  £,  is  the  energy  of  molecule  i,  ej  (denot¬ 
ed  as  an  “atom  equivalent”)  is  defined  as  ej  = 
(Ej  —  Xj),  Ej  is  the  energy  of  an  atom  j  that  is 
a  component  of  molecule  i,  rij  is  the  number  of 
j  atoms  in  molecule  i,  andxy  is  the  correction  for 
atom  j  at  the  level  of  theory  used.  The  atom 
equivalents  are  determined  through  least- 
squares  fitting  of  Eq.  3  to  experimental  heats  of 
formation  and  quantum  mechanical  energies 
for  several  representative  molecules.  The  atom 
equivalents  correct  for  the  error  inherent  in  the 
calculation  and  in  some  cases,  for  the  tempera¬ 
ture  and  zero-point  energy.  Atom  equivalent 
schemes  are  popular  alternative  methods  for 
predicting  heats  of  formation  since  they  do  not 
require  high-level  treatment  of  electron  corre¬ 
lation  or  experimental  input  once  the  atom 
equivalents  have  been  determined  [16-20].  This 
approach  was  recently  applied  to  small  hydro¬ 
carbons  [16]  and  a  variety  of  organic  com¬ 
pounds  [20]  using  DFT.  The  Mole  et  al.  study 
found  that  B3LYP/6-31G*  predictions  of  heats 
of  formation  of  23  hydrocarbons  had  a  root 
mean  square  (rms)  deviation  of  1.72  kcal/mol, 
with  the  maximum  difference  between  experi¬ 
ment  and  prediction  being  6.2  kcal/mol  [16]. 
Heats  of  formation  calculated  using  the 
6-311 +G**  basis  set  improved  the  accuracy  of 
the  predictions  to  a  rms  deviation  from  experi¬ 
ment  of  1  kcal/mol,  and  a  maximum  deviation  of 
2.34  kcal/mol.  Habibollahzadeh  et  al.  [20]  used 
a  similar  procedure  to  predict  the  gas-phase 
heats  of  formation  of  54  organic  compounds 
using  the  Becke  exchange  and  Perdew  correla¬ 
tion  functionals  [21,  22]  and  the  6-31G(d,p) 
basis  set  [14].  The  results  have  an  average 
deviation  from  experiment  of  3  kcal/mol. 

We  will  follow  the  approach  given  in  Ref.  16 
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to  determine  atom  equivalents  of  carbon,  nitro¬ 
gen,  hydrogen,  and  oxygen,  the  atomic  species 
in  most  energetic  materials.  To  do  this,  we  first 
calculate  a  set  of  energies  that  correspond  to 
optimized  structures  of  energetic  molecules  for 
which  heats  of  formation  have  been  measured. 
We  then  determine  the  atom  equivalents 
through  least-squares  fitting  of  Eq.  3  using  the 
set  of  quantum  mechanical  energies  and  mea¬ 
sured  heats  of  formation.  The  ideal  approach  is 
to  determine  atom  equivalents  that  do  not 
depend  on  bond  order  or  group  component. 
However,  there  are  several  energetic  molecules 
distinguished  by  different  functional  groups  and 
different  bond  orders,  but  that  have  the  same 
molecular  formula.  Because  of  this,  we  have 
determined  the  atom  equivalents  for  atoms 
involved  in  a  single  Lewis  structure  (two-elec¬ 
tron  bond)  or  a  multiple  Lewis  structure  (great¬ 
er  than  a  two-electron  bond,  such  as  seen  in 
aromatic  groups  and  nitro  groups).  The  heat  of 
formation  as  given  in  Eq.  3  is  therefore  deter¬ 
mined  through  seven  atom  equivalents,  four 
representing  atoms  involved  in  single  bonds 
(denoted  as  C,  H,  N,  and  O),  and  three  repre¬ 
senting  atoms  in  multiple  bonds  (denoted  as  C', 
N',  O').  The  atom  equivalents  were  determined 
by  least-squares  fitting  of  experimental  heats  of 
formation  and  B3LYP/6-31G*  molecular  ener¬ 
gies  for  35  molecules  to  Eq.  3. 

Although  Eq.  3  has  been  shown  to  be  very 
accurate  in  the  prediction  of  gas-phase  heats  of 
formation,  often  the  standard  state  of  the  ma¬ 
terial  of  interest  corresponds  to  the  condensed 
phase.  Thus,  it  is  important  to  be  able  to  predict 
the  condensed-phase  heats  of  formation.  Con¬ 
densed-phase  heats  of  formation  can  be  deter¬ 
mined  using  the  gas-phase  heat  of  formation 
and  heat  of  phase  transition  (either  sublimation 
or  vaporization)  according  to  Hess’  law  of  con¬ 
stant  heat  summation  [11]: 

AH(Solid)  =  A//(Gas) 

-  AH(Sublimation)  (4) 

AH(Liquid)  =  A//(Gas) 

—  AFr(  Vaporization).  (5) 

To  evaluate  Eqs.  4  and  5  using  theoretical 
predictions  alone,  tools  must  be  developed  to 


predict  enthalpies  of  phase  change.  Such  tools 
have  been  developed  by  Politzer  and  coworkers 
for  different  types  of  organic  compounds  with  a 
significant  degree  of  success  [23-29]. 

In  a  series  of  studies  that  investigate  the 
relation  of  bulk  properties  of  a  niaterial  with 
molecular  properties  [23-29],  Politzer  and  co¬ 
workers  have  clearly  established  that  correla¬ 
tions  exist  between  the  electrostatic  potential  of 
a  molecule  and  its  condensed-phase  properties, 
including  the  heats  of  sublimation  [27]  and 
vaporization  [24,  26].  Functional  relationships 
based  on  such  correlations  can  be  determined 
through  the  following  procedure:  First,  Politzer 
et  al.  recommend  determining  the  low-energy 
structure  of  a  molecule  through  quantum  me¬ 
chanical  prediction,  then  calculating  the  molec¬ 
ular  surface  area  (SA)  for  this  structure  [23-29]. 
In  this  approach,  the  SA  is  defined  to  be  that 
corresponding  to  the  0.001  electrons/bohr^  iso¬ 
surface  of  the  electron  density.  Next,  the  elec¬ 
trostatic  potential  for  this  isosurface  is  used  to 
generate  two  statistically  based  quantities, 
and  V.  Oj-o,  is  described  as  an  indicator  of  the 
variability  of  the  electrostatic  potential  on  the 
molecular  surface,  and  v  is  interpreted  as  show¬ 
ing  the  degree  of  balance  between  the  positive 
and  negative  potentials  on  the  molecular  sur¬ 
face.  From  these  quantities,  the  heat  of  vapor¬ 
ization  can  be  represented  as  [24,  26]: 

A//(  Vaporization)  =  a  V(SA)  +  byjo^^  +  c 

(6) 

where  a,  b,  and  c  are  fitting  parameters.  Simi¬ 
larly,  the  heat  of  sublimation  can  be  represented 
as  [27]: 

A//(Sublimation)  =  a(SA)^  +  b  yjo'Toi'^  +  ^ 

(7) 

The  parameters  a,  b,  and  c  in  Eqs.  6  and  7  are 
determined  from  least-squares  fitting  to  reliable 
values  of  the  enthalpies  of  phase  change. 
Politzer  et  al.  applied  this  procedure  to  develop 
a  predictive  tool  for  the  heats  of  sublimation  of 
34  organic  compounds  [27].  The  predictions 
using  the  tool  had  a  standard  deviation  of  2.5 
kcal/mol  from  experimental  values.  Politzer  et 
al.  also  used  this  tool  with  one  developed  to 
calculate  gas-phase  heats  of  formation  to  pre- 


448 


B.  M.  RICE  ET  AL. 


diet  solid-phase  heats  of  formation  for  five 
compounds  [27].  In  this  application,  the  average 
deviation  from  experiment  was  2.8  kcal/mol.  In 
a  similar  application,  Politzer  and  Murray  pre¬ 
dicted  heats  of  vaporization  using  Eq.  5  for  41 
compounds  that  have  a  standard  deviation  of 
0.6  kcal/mol  [26]. 

We  have  followed  the  Politzer  et  al.  approach 
to  develop  tools  to  predict  the  heats  of  sublima¬ 
tion  and  vaporization  of  energetic  materials, 
then  combined  these  with  Eq.  3  to  predict  solid 
and  liquid  heats  of  formation.  The  only  signifi¬ 
cant  difference  in  the  work  presented  here  and 
that  presented  by  Politzer  and  coworkers  is  the 
level  of  theory  and  the  set  of  molecules  used  in 
parameterizing  Eqs.  6  and  7. 

COMPUTATIONAL  DETAILS 

Generalized  gradient  approximation  (GGA) 
DFT  [4-7]  geometry  optimizations  of  all  species 
reported  herein  were  performed  using  the 
6-31G*  basis  set  [14]  and  the  hybrid  B3LYP  [9, 
10]  density  functional.  This  modest  level  of 
theory  was  chosen  due  to  the  size  of  some  of  the 
molecules  used  in  the  study.  The  6-31G*  basis 
set  has  been  shown  to  be  reasonably  accurate 
when  used  with  the  B3LYP  density  functional 
[16,  30].  The  B3LYP  density  functional  has 
been  shown  to  reproduce  experimental  proper¬ 
ties  and  is  commonly  thought  to  be  one  of  the 
most  reliable  of  the  available  density  functionals 
[1,  8,  31].  The  calculations  were  performed 
using  the  Gaussian  94  (G94)  suite  of  quantum 
chemistry  programs  [32].  All  calculations  were 
subject  to  the  default  settings  of  G94. 

The  atom  equivalents  in  Eq.  3  are  determined 
by  fitting  to  experimental  gas-phase  heats  of 
formation  and  corresponding  quantum  mechan¬ 
ical  energies  of  optimized  structures  of  several 
energetic  molecules.  The  quantum  mechanical 
energies  of  the  molecules  whose  names  are 
followed  by  a  “g”  in  Table  1  were  used  in  the 
fitting  of  Eq.  3.  Many  of  the  molecules  used  in 
the  fitting  are  large  polyatomic  molecules,  and 
most  likely  have  several  different  stable  confor¬ 
mations  that  are  similar  in  energy.  We  did  not 
perform  a  search  in  conformation  space  for 
each  of  these  molecules  to  determine  the  mo¬ 
lecular  structure  associated  with  the  global  min¬ 


imum,  but  have  assumed  that  the  energy  of  the 
local  minima  are  within  a  few  kcal/mol  of  the 
global  minimum.  Eqs.  6  and  7  are  also  param¬ 
eterized  through  least-squares  fitting  to  experi¬ 
mental  data  and  information  from  the  quantum- 
mechanicaUy-derived  electrostatic  potential  of 
individual  molecules,  as  described  previously. 
Equations  6  and  7  were  fitted  using  27  and  36 
sets  of  experimental  data  for  heats  of  vaporiza¬ 
tion  and  sublimation,  respectively.  The  mole¬ 
cules  used  in  parameterizing  Eqs.  6  and  7  are 
also  given  in  Table  1.  Their  names  are  followed 
by  a  “v”  or  “s”  to  denote  that  they  were  used  for 
parameterizing  the  equations  for  the  heat  of 
vaporization  or  sublimation,  respectively. 

RESULTS  AND  DISCUSSION 

Atom  equivalents  obtained  through  least- 
squares  fitting  of  Eq.  3  are  given  in  Table  2.  A 
comparison  of  the  experimental  and  the  pre¬ 
dicted  gas-phase  heats  of  formation  using  Eq.  3 
is  given  in  Table  3.  A  visual  comparison  of  the 
predictions  with  experiment  is  provided  in  Fig. 
1.  The  rms  deviation  of  the  predictions  from 
experiment  is  3.1  kcal/mol.  The  maximum  devi¬ 
ation  of  the  predictions  from  the  experimental 
values  is  7.3  kcal/mol  for  azidomethylbenzene. 
The  next  largest  difference  occurs  for  cyanogen 
azide  (6.2  kcal/mol).  This  is  not  surprising,  since 
only  two  azido-  compounds  were  used  in  the 
fitting. 

The  best-fit  parameters  for  Eqs.  6  and  7  are 
given  in  Table  2,  and  the  predicted  heats  of 
vaporization  and  sublimation  are  given  in  Table 
3  for  comparison  with  the  experimental  infor¬ 
mation  used  in  the  fitting.  Figures  2  and  3  show 
the  results  of  the  fits  for  the  heats  of  vaporiza¬ 
tion  and  sublimation,  respectively.  The  rms  de¬ 
viation  of  the  predictions  from  the  experimental 
heats  of  vaporization  is  1.7  kcal/mol.  The  max¬ 
imum  deviation  from  experiment  is  6.1  kcal/mol 
for  nitroglycerin,  with  the  next  largest  differ¬ 
ence  being  2.9  kcal/mol  for  trinitromethane. 
The  rms  deviation  of  the  predictions  from  ex¬ 
perimental  heats  of  sublimation  is  3.6  kcal/mol, 
with  the  maximum  deviation  being  12.4  kcal/mol 
for  hexanitroethane.  The  next  largest  difference 
between  experiment  and  prediction  is  5.3  kcal/ 
mol  for  dinitromethylbenzene. 
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TABLE  1 


B3LYP/6-31G*  Molecular  Properties  Used  in  Eqs.  3,  6,  and  7 


Name 

Expt." 

c 

H 

N 

o 

N' 

O' 

Absolute  Energy 
(hartrees) 

SA  (A^) 

V 

oto,  (kcal/ 
mol)2 

Cyanogen  azide 

g 

0 

0 

0 

0 

1 

4 

0 

-257.002988577 

97.646679 

0.249860 

207.859558 

Tetranitromethane 

g.  V 

1 

0 

0 

0 

0 

4 

8 

-858.413321934 

166.035026 

0.052589 

82.899727 

Trinitromethane 

V,  s 

1 

1 

0 

0 

0 

3 

6 

-653.956381558 

144.477765 

0.043476 

259.188354 

Hexanitroethane 

g.  s 

2 

0 

0 

0 

0 

6 

12 

-1306.656763240 

212.868693 

0.064430 

119.363091 

Dinitromethane 

g.  V 

1 

2 

0 

0 

0 

2 

4 

-449.487256631 

118.551923 

0.100087 

240.648117  . 

Methyl  nitrite 

g»  V 

1 

3 

0 

1 

0 

1 

1 

-245.007962059 

89.882874 

0.249936 

56.139923 

Nitromethane 

g»  V 

1 

3 

0 

0 

0 

1 

2 

-245.009330800 

89.818926 

0.246111 

108.135544 

Methyl  nitrate 

g^v 

1 

3 

0 

1 

0 

1 

2 

-320.189434600 

100.572911 

0.249294 

79.381935 

Nitroguanidine 

0 

4 

2 

0 

1 

2 

2 

-409.855059996 

124.063733 

0.230870 

435.239868 

Nitroethane 

V 

2 

5 

0 

0 

0 

1 

2 

-284.328086959 

110.468975 

0.229348 

103.713966 

Ethyl  nitrite 

2 

5 

0 

1 

0 

1 

1 

-284.325006418 

115.069420 

0.213790 

60.119099 

Ethyl  nitrate 

g’ V 

2 

5 

0 

1 

0 

1 

2 

-359.508395169 

121.367742 

0.241344 

69.094345 

Dimethylnitramine 

g.s 

2 

6 

1 

0 

0 

1 

2 

-339.656549984 

123.494464 

0.200675 

140.171310 

Nitroglycerin 

g,  V 

3 

5 

0 

3 

0 

3 

6 

-958.168013834 

218.847052 

0.111295 

132,965546 

TTT  (Hexahydro-1,3,5- 

g»s 

3 

6 

3 

0 

0 

3 

3 

-671.860924393 

181.137911 

0.211136 

162.981583 

trinitroso-l,3,5-triazine) 

RDX  (Hexahydro-1,3,5- 

g^S 

3 

6 

3 

0 

0 

3 

6 

-897.408901632 

201.594799 

0.149852 

196.611298 

trinitrotriazine) 

l-Nitropropane 

V 

3 

7 

0 

0 

0 

1 

2 

-323.641585437 

132.373802 

0.199392 

92.913933 

2-Nitropropane 

V 

3 

7 

0 

0 

0 

1 

2 

-323.645738999 

129.743605 

0.182169 

91.661804 

Propyl  nitrite 

g»  V 

3 

7 

0 

1 

0 

1 

1 

-323.638914025 

136.474131 

0.181451 

56.687702 

n-Propyl  nitrate 

V 

3 

7 

0 

1 

0 

1 

2 

-398.820582169 

140.472688 

0.238864 

70.809593 

3,4-Furazandimethanol 

g.  V 

2 

4 

0 

3 

2 

4 

4 

-900.030164899 

211.827748 

0.127652 

123.151146 

dinitrate 

1,4-Dinitrosopiperazine 

g.  S 

4 

8 

2 

0 

0 

2 

2 

-526.543460553 

168.541711 

0.225159 

111.991249 

1,4-Dinitropiperazine 

g.  S 

4 

8 

2 

0 

0 

2 

4 

-676.910151215 

184.191543 

0.149480 

110.465118 

HMX  (Octahydro-1,3^,7- 

4 

8 

4 

0 

0 

4 

8 

-1196.545948910 

245.151259 

0.188715 

196.189636 

tetranitro-1,3,5,7- 

tetrazocine) 

2-Methyl-2-nitropropane 

g 

4 

9 

0 

0 

0 

1 

2 

-362.961528577 

144.467749 

0.159268 

97.889603 

1-Nitrobutane 

V 

4 

9 

0 

0 

0 

1 

2 

-362.954119935 

151.853360 

0.180655 

91.500481 

2-Nitrobutane 

V 

4 

9 

0 

0 

0 

1 

2 

-362.958529381 

148.815487 

0.145739 

86.887459 

t-Butyl  nitrite 

g»  V 

4 

9 

0 

1 

0 

1 

1 

-362.958313184 

149.240879 

0.161770 

66.320312 

n-Butyl  nitrite 

g»  V 

4 

9 

0 

1 

0 

1 

1 

-362.952557150 

157.775266 

0.178047 

57.120720 

PETN  (Tetranitrate 

s 

5 

8 

0 

4 

0 

4 

8 

-1316.468455720 

291.516975 

0.104211 

116.481842 

pentaerythritol) 

1-Nitropiperidine 

g 

5 

10 

1 

0 

0 

1 

2 

-456.398610158 

163.824522 

0.095410 

123.016701 

1,3,5-Trinitrobenzene 

s 

0 

3 

0 

0 

6 

3 

6 

-845.736696893 

202.532990 

0.188445 

108.958405 

2,4,6-T  rinitrophenol 

0 

3 

0 

1 

6 

3 

6 

-920.954678338 

207.730401 

0.198037 

135.451736 

2,4,6-Trinitroresorcinol 

0 

3 

0 

2 

6 

3 

6 

-996.177452790 

214.610239 

0.194153 

150.847977 

l-Nitro-2-nitrosobenzene 

s 

0 

4 

0 

0 

6 

2 

3 

-565.992564774 

164.664846 

0.249159 

103.950478 

l-Nitro-3-nitrosobenzene 

s 

0 

4 

0 

0 

6 

2 

3 

-566.036632071 

168.137290 

0.249814 

99.408859 

1,2-Dinitrobenzene 

0 

4 

0 

0 

6 

2 

4 

-641.229491455 

175.485624 

0.247027 

125.204185 

1,3-Dinitrobenzene 

s 

0 

4 

0 

0 

6 

2 

4 

-641.246252172 

176.131173 

0.248380 

97.582031 

1,4-Dinitrobenzene 

0 

4 

0 

0 

6 

2 

4 

-641.246008400 

175.928160 

0.249999 

84.098793 

2,4-Dinitrophenol 

0 

4 

0 

1 

6 

2 

4 

-716.473897710 

182.360129 

0.249977 

103.774368 

2,6-Dinitrophenol 

0 

4 

0 

1 

6 

2 

4 

-716.461450000 

180.503403 

0.246903 

124,667587 

2,4,6-Trinitroaniline 

0 

4 

1 

0 

6 

3 

6 

-901.105946802 

209.342650 

0.218842 

125.144867 

Nitrosobenzene 

g 

0 

5 

0 

0 

6 

1 

1 

-361.539779255 

142.410932 

0.199211 

146.791779 

Nitrobenzene 

g>  V 

0 

5 

0 

0 

6 

1 

2 

-436.750579273 

150.299993 

0.215926 

154.002136 

2-Nitrophenol 

g^  S 

0 

5 

0 

1 

6 

1 

2 

-511.976038167 

156.500748 

0.236806 

114.630844 

3-Nitrophenol 

g.  S 

0 

5 

0 

1 

6 

1 

2 

-511.966822213 

159.153839 

0.232478 

248.392502 

4-Nitrophenol 

g’ s 

0 

5 

0 

1 

6 

1 

2 

-511.969008406 

159.159240 

0.234619 

305.803284 

m-Nitroaniline 

g.  S 

0 

6 

1 

0 

6 

1 

2 

-492.104971819 

164.287003 

0.248321 

256.141357 
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TABLE  1 


(Continued) 


Name 

Expt.'' 

c 

H 

N 

0 

C' 

N' 

O' 

Absolute  Energy 
(hartrees) 

SA  (A^) 

V 

(kcal/ 

mol)2 

o-Nitroaniline 

s 

0 

6 

1 

0 

6 

1 

2 

-492.109521689 

160.796544 

0.247430 

232.851974 

TATB  (2,4,6-trinitro-l,3,5- 

0 

6 

3 

0 

6 

3 

6 

-1011.833118600 

221.751885 

0.249990 

116.236404 

benzenetriamine) 

DNPN  [N-Nitrobis-(2,2- 

g 

6 

10 

1 

0 

0 

5 

10 

-1314.872258970 

276.518965 

0.206557 

141.827118 

dinitropropyl)amine] 

TNT  (trinitrotoluene) 

g.s 

1 

5 

0 

0 

6 

3 

6 

-885.045499312 

217.011451 

0.215058 

94.741478 

2-Methoxy-l,3,5- 

s 

1 

5 

0 

1 

6 

3 

6 

-960.240793345 

226.444172 

0.221238 

96.778053 

trinitrobenzene 

Tetryl  (N-methyl-N,2,4,6- 

1 

5 

1 

0 

6 

4 

8 

-1144.856680870 

251.412792 

0.178351 

141.468536 

tetranitroaniline) 

l-Methyl-2,4- 

g.s 

1 

6 

0 

0 

6 

2 

4 

-680.561775882 

192.841213 

0.246866 

90.642029 

dinitrobenzene 

2-Methyl-l,3- 

1 

6 

0 

0 

6 

2 

4 

-680.553580234 

190.525356 

0.248539 

98.388885 

dinitrobenzene 

Dinitromethylbenzene 

g.s 

1 

6 

0 

0 

6 

2 

4 

-680.543248757 

196.337841 

0.237477 

122.123009 

2-Methyl-4,6-dinitrophenol 

1 

6 

0 

1 

6 

2 

4 

-755.792238817 

201.448137 

0.236338 

91.450188 

Azidomethylbenzene 

g.v 

1 

7 

0 

0 

6 

3 

0 

-435.144931168 

179.499171 

0.232174 

65.953873 

Nitromethylbenzene 

g.v 

1 

7 

0 

0 

6 

1 

2 

-476.061676752 

172.640031 

0.199769 

129.356613 

l-Methyl-4-nitrobenzene 

g.s 

1 

7 

0 

0 

6 

1 

2 

-476,069640689 

171,147817 

0.171034 

158.963928 

2,4,6-Trinitrometaxylene 

2 

7 

0 

0 

6 

3 

6 

-924.357319040 

232.435743 

0.244814 

73.367279 

l,3-Dimethyl-2- 

g.v 

2 

9 

0 

0 

6 

1 

2 

-515.379313077 

185.404942 

0.162293 

130.591660 

nitrobenzene 

l,3,5-Trimethyl-2,4,6- 

3 

9 

0 

0 

6 

3 

6 

-963.668410824 

248.796721 

0.248856 

57,152065 

trinitrobenzene 

1-Nitronaphthalene 

s 

0 

7 

0 

0 

10 

1 

2 

-590.387677028 

194.108673 

0.198561 

159.976227 

2,2',4,4',6,6'- 

g 

0 

6 

0 

0 

14 

6 

12 

-1767.644192350 

381.660470 

0.166741 

119.739990 

Hexanitrostilbene 

e-Hexaazaisowurtzitane 

6 

6 

6 

0 

0 

6 

12 

-1791.180675600 

312.693403 

0.065415 

230.603210 

/3-Hexaazaisowurtzitane 

6 

6 

6 

0 

0 

6 

12 

-1791.183184600 

314.537686 

0.071639 

213.575226 

Dinitrate  diethylene  glycol 

4 

8 

0 

3 

0 

2 

4 

-793.010912406 

218.336681 

0.243205 

69.461563 

Symbols  denote  experimental  values  are  available  for  gas-phase  heats  of  formation  (g),  heats  of  vaporization  (v),  and  heats 
of  sublimation  (s)  and  are  used  in  parameterizing  Eqs.  3,  6,  and  7. 


Table  3  also  contains  comparisons  of  predic¬ 
tions  and  experimental  values  for  liquid  and 
solid  heats  of  formation.  There  are  41  experi¬ 
mental  values  of  liquid  heats  of  formation  that 


correspond  to  24  molecules.  These  are  com¬ 
pared  with  predictions  using  Eqs.  3  and  6. 
Figure  4  provides  a  visual  comparison  of  the 
predictions  and  experiment.  The  rms  deviation 


TABLE  2 

Atom  Equivalents  and  Parameters  for  Eqs.  6  and  7 


Equation  3 


Atom  Equivalents  e  (hartrees)  Equation  6  Equation  7 


C 

-38.121621 

a  (kcal/mol-A"^) 

1.818689 

a  (kcal/mol- A""^) 

4.234303  X  10 

H 

-0.592039 

b  (kcal/mol) 

1,3321583 

b  (kcal/mol) 

2.5793785 

N 

-54.774096 

c  (kcal/mol) 

-16.142460 

c  (kcal/mol) 

-6.7335407 

O  -75.161771 

a  -38.121380 

N'  -54.765886 

O'  -75.157348 


Heats  of  Formation,  Heat  of  Vaporization,  and  Heat  of  Sublimation  (kcal/mol)' 
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EXPERIMENTAL  (kcal/mol) 


Fig.  1.  Calculated  gas-phase  heats  of  formation  versus 
experimental  values  for  35  energetic  molecules.  The  solid 
line  represents  exact  agreement  between  the  predictions 
and  experiment. 

of  the  predicted  liquid  heats  of  formation  from 
experiment  is  3.3  kcal/mol,  and  the  maximum 
deviation  from  experiment  is  9.3  kcal/mol  for 
azidomethylbenzene,  followed  by  one  of  the 
values  reported  for  trinitromethane.  Note  that 
three  values  are  reported  for  trinitromethane; 
these  range  from  —7.68  to  —18.63  kcal/mol.  The 
predicted  value  is  —10.8  kcal/mol.  The  only 
molecular  species  in  Fig.  4  for  which  no  exper¬ 
imental  information  was  used  in  parameterizing 
Eqs.  3  or  6  is  dinitrate  diethylene  glycol;  the 
predicted  liquid  heat  of  formation  is  within  3 
kcal/mol  of  the  experimental  result. 


Fig.  2.  Same  as  Fig.  1,  except  for  heats  of  vaporization. 
Comparisons  are  made  using  27  experimental  values. 


Fig.  3.  Same  as  Fig.  1,  except  for  heats  of  sublimation. 
Comparisons  are  made  using  36  experimental  values. 


Seventy-five  experimental  values  for  solid 
heats  of  formation  corresponding  to  44  mole¬ 
cules  are  given  in  Table  3  for  comparison  with 
predictions  using  Eqs.  3  and  7.  Figure  5  provides 
a  visuzil  comparison  between  experiment  and 
the  predictions.  The  rms  deviation  of  the  pre¬ 
dictions  from  experiment  is  9.0  kcal/mol,  and 
the  maximum  deviation  of  35.4  kcal/mol  corre¬ 
sponds  to  one  of  the  reported  values  for  trini- 
troresorcinol.  As  evident  from  Table  3,  there 


EXPERIMENTAL  (kcal/mol) 


Fig.  4.  Same  as  Fig.  1,  except  for  liquid-phase  heats  of 
formation.  Comparisons  are  made  using  41  experimental 
values. 


PREDICTING  HEATS  OF  FORMATION 


Fig.  5.  Same  as  Fig.  1,  except  for  solid-phase  heats  of 
formation.  Comparisons  are  made  using  75  experimental 
values. 


are  several  species  for  which  notably  different 
values  are  reported  for  the  heat  of  formation. 
The  values  reported  for  solid  heats  of  formation 
for  trinitroresorcinol,  for  example,  range  from 
—103.5  to  —129.8  kcal/mol.  The  predicted  heat 
of  formation  deviates  from  one  of  the  reported 
values  by  35.4  kcal/mol,  which  is  the  maximum 
deviation  of  all  of  the  predictions.  The  same 
prediction,  however,  is  within  9  kcal/mol  of  one 
of  the  other  reported  values  for  trinitroresor¬ 
cinol.  The  next  largest  deviation  of  the  predic¬ 
tions  from  experiment  is  for  hexanitrostilbene, 
which  is  —26  kcal/mol. 

Eighteen  of  these  molecular  species  were  not 
used  in  the  parameterization  of  Eqs.  3  or  7,  and 
thus  the  results  for  these  can  be  considered  as 
completely  predicted,  rather  than  reflecting  the 
degree  of  the  goodness  of  the  fitting  of  Eqs.  3 
and  7.  Figure  6  shows  the  comparison  of  the 
predictions  from  experiment.  The  rms  deviation 
of  the  predictions  from  experiment  for  these  18 
molecular  species  is  10.3  kcal/mol,  respectively, 
and  the  maximum  deviation  from  experiment  is 
11.5  kcal/mol  for  Tetryl.  Four  of  the  18  species 
have  multiple  values  reported  for  heats  of  for¬ 
mation  that  differ  from  one  another  by  more 
than  5  kcal/mol.  The  values  reported  for  HMX 
differ  by  6.5  kcal/mol;  those  of  2,4,6-trinitrores- 
orcinol  differ  by  26.3  kcal/mol;  those  of  2,4,6 
trinitroaniline  differ  by  10.3  kcal/mol;  and  those 
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EXPERIMENTAL  (kcal/mol) 

Fig.  6.  Same  as  Fig.  5,  except  it  shows  a  comparison 
between  predictions  and  experiment  for  18  molecular  spe¬ 
cies  from  which  no  information  was  used  in  developing  the 
predictive  tools  described  in  this  work. 

of  TATB  differ  by  19.1  kcal/mol.  Such  dispari¬ 
ties  in  experimental  measurement  contribute  to 
the  collective  rms  deviation  of  the  predictions 
from  experiment. 

CONCLUSIONS 

A  series  of  computational  tools  have  been  de¬ 
veloped  that  use  only  quantum  mechanical  in¬ 
formation  to  predict  gas-  and  condensed-phase 
heats  of  formation,  heats  of  vaporization,  and 
heats  of  sublimation.  Quantum  mechanical  en¬ 
ergies  of  molecules  are  converted  to  gas-phase 
heats  of  formation  using  the  method  of  atom 
equivalents  [16].  The  predictions  for  35  mole¬ 
cules  have  a  rms  deviation  from  experimental 
results  of  3.1  kcal/mol.  Surface  electrostatic 
potentials  of  individual  molecules  were  gener¬ 
ated  and  used  in  empirical  equations  for  heats 
of  sublimation  and  vaporization  as  recom¬ 
mended  by  Politzer  and  coworkers  [23-29]. 
These  functions  were  parameterized  using  27 
and  36  sets  of  experimental  values  for  the  heats 
of  vaporization  and  sublimation,  respectively. 
The  rms  deviation  of  the  heats  of  vaporization 
calculated  using  this  function  is  1.7  kcal/mol. 
The  rms  deviation  of  heats  of  sublimation  from 
experimental  values  is  3.6  kcal/mol.  The  equa- 
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tions  for  heats  of  vaporization  and  sublimation 
are  used  with  the  atom-equivalent  method  for 
calculating  the  gas-phase  heats  of  formation  to 
predict  liquid  and  solid  heats  of  formation.  The 
liquid  heats  of  formation  have  a  rms  deviation 
from  41  experimental  values  of  3.3  kcal/mol. 
The  rms  deviation  of  the  solid  heats  of  forma¬ 
tion  from  75  experimental  values  is  9.0  kcal/mol. 
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